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Abstract 

We review some recent applications of QCD sum rules to B- 
decays. These include the determination of the leptonic con- 
stant Jb and of the semileptonic transition amplitudes of the 
S-meson into negative and positive parity charmed states. 
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1. Introduction 



Considerable interest has been given in the past years to applications of 
QCD sum rules to B (and D) physics, resulting in a large number of estimates 
for the B and D transition matrix elements. 

The formalism consists in relating via dispersion relations low energy pa- 
rameters, such as hadronic masses and couplings, to the QCD dominated short- 
distance Operator Product Expansion (OPE) of current correlators [1,2]. Such 
an expansion is assumed to hold in the presence of non-perturbative QCD ef- 
fects, which are parametrized by a set of fundamental vacuum condensates, 
describing the breaking of asymptotic freedom, and inducing power corrections 
in the correlators. These condensates are fundamental constants which can- 
not be computed from perturbation theory: however, they are "universal" and, 
once fitted from channels where the physics is known, their values can be used 
to make predictions for the other channels. Thus, being fully field-theoretic and 
relativistic, and incorporating the fundamental features of QCD, this method 
aims to be a means to interface from first principles hadronic observables to 
the basic quark and gluon QCD Lagrangian. From this point of view, it should 
represent a convenient framework to study the properties of heavy-light quark 
systems and of quarkonia, which avoids the notion of constituent quark bound 
state wavefunction and should thus reduce the model dependence to a minimum. 

Most recently, the attention has been focussed on the applications of QCD 
sum rules to the semileptonic and the purely leptonic B- and D-meson decays, 
motivated on one side by the observed features of the data [3,4] (also in the 
perspectives of increased precision and of rare decays measurements, as allowed 
by planned charm and beauty factories), and on the other side by the theoretical 
achievements offered in principle by the heavy quark effective theory (HQET) 
[5]. In this regard, it turns out that QCD sum rules allow to study the behaviour 
of amplitudes for qQ mesons for wlq — > oo (Q represents the heavy quark, q the 
light one). Therefore they can be used to verify the predictions of spin-flavor 
heavy quark symmetry, and to test quantitatively the resulting scaling laws 
relating the different amplitudes as well as the validity of approximations based 
on the 1/rriQ expansions. Accordingly, the sequel will be devoted to applications 
of QCD sum rules to several aspects of B physics. Specifically, we shall briefly 
review the predictions for the leptonic constant analyzing in particular the 
determination of the beauty quark mass m&, and of the semileptonic S-decays 
to negative and to positive parity charmed mesons. 
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2. The leptonic decay constant fs 

Leptonic decay constants, governing the decays D — > /jo;, B — > etc., are 
defined as the vacuum-to-meson matrix elements: 

< 0|Q 7M 7 5 g|P(p) >= if P (1) 

with Q = c,b the heavy quarks and q = u,d, s the light ones, so that P = B, D. 
The theoretical interest stems from their significant role as normalizations of 
numerous predictions for heavy meson transitions. In particular, by dominating 
the description of B — B mixing, f B strongly affects the determination of the 
CKM matrix, with important implications on CP violation in the beauty sector 
[6] . Furthermore, of great interest is the scaling behaviour for very large (infinite) 

const 

JP = , > \ l ) 

modulo logs, which reflects the connection with the constituent quark wavefunc- 
tion, but in fact results from the heavy quark symmetry and accordingly is a 
general consequence of QCD. 

The present indications are f D < 310 MeV [7] and f Da = (232±45±20±48) 
MeV [8], whereas no information is available on /#. To get a feeling on the 
sensitivity required to measure /#, for (V^/Kbl = 0.1 with \V c b\ = 0.045 and 
f B = 1.5/; (fn = 132 MeV) one finds BR(B -> fiu) = 5.1 x 10" 7 and BR(B -> 

TV) = 1.1 X 10" 4 . 

Let us now briefly describe the method of QCD sum rules as applied to /#. 
One starts from the two-point correlator at euclidean Q 2 = —q 2 : 

U 5 (q 2 ) = i J d 4 xexp(iqx) < 0|T (J 5 (x)J 5 (0) f ) |0 >, (3) 

where J$(x) =: Q{x)i^q{x) : is the quark bilinear operator which interpolates 
the pseudoscalar meson field, and one can neglect the light quark mass m q 
with respect to tuq. For other mesonic states, the appropriate interpolating 
fields must be used in (3), e.g. V M =: Q^^q : for the vector particle B* , etc. 
In QCD H 5 (q 2 ) satisfies a twice-subtracted dispersion relation. Actually, to 
improve the convergence, emphasize the contribution of the hadronic ground 
state one is interested in, and get rid of unknown subtraction constants, suitable 
weight functions are used in the dispersive integrals. One is thus led to consider 
"moments" , each corresponding to a different version of QCD sum rules, such 
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as the Hilbert moments at Q 2 = 0, leading, for n = 1, 2 . . ., to the finite energy 
sum rules (FESR): 

(-l) n / d \ n+1 r°° li 

M » (0) s fcm W) n ^ >Uo = / ^^M; (4) 

and the exponential moments 

Z" 00 1 
115(a) = / dsex.p(—as)—lmll 5 (s), (5) 
Jo ' ^~ 

where 1 / \fa is a mass scale. A review of other versions can be found in [2] . 

The hadronic spectral function on the right sides of (4) and (5) is usually 
parametrized in terms of the ground state meson pole plus a "continuum" of 
higher states, starting at some threshold sq and modelled by the asymptotic 
freedom (AF) expression: 

-ImU 5 (s)\ HAD = &^8(s - m 2 P ) + -ImU 5 (s) \ 6(s - s ). (6) 

7V 712 q TV 



In (6), to two loops [9]: 

1 r^rr ,„M _ 3 ,....2 - •<•>- ! -, , fas 

3tt 



1 T „ / M 3 2 (1-x) 2 r 



18 nT . f -x 
— - 2Lz 2 



1 - x 



- In I — ^— 1 In I — — I + I - - zr— x ) In 



1 — x J \1 — x J \2 1 — x J \1 — x r 

(7) 



1 . / 1 
H In 



1 — x \ 1 — x 
where x = ijiq/s. 

The left sides of (4) and (5) can be evaluated by means of the OPE, in 
terms of short-distance perturbative QCD, determined by the integral of (7) 
from m,Q up to infinity, plus non-perturbative power corrections, determined by 
quark and gluon operator vacuum condensates, ordered for increasing dimen- 
sion. These vacuum expectation values essentially allow the extrapolation of 
the asymptotic freedom contribution, completely known in terms of a s and uiq, 
down to "moderate" mass scales, relevant to hadrons. 

Limiting to operators of dimensionality up to six, combining (4), (6) and 
(7) one obtains the FESR, in the form (m q — > 0): 

fp 2 f ds 1 



/ds 1 
g n+2 ~^ m ^ 5 I AF 
+ M n (0)\ Np , (8) 
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where the non-perturbative part is given by: 



M n (0)| 



-m Q (qq) 



NP 



m 



2n+2 

Q 



1 " To - \ - 7 ( w + 2 )(n + l)-f 
12irm Q (qq) 4 m 2 



Q 



~(w + 2)(n 2 + lOra + 9)7ra s/ o^ 
81 



(9) 



Alternatively, combining (5), (6) and (7) one obtains the exponential (also called 
Laplace or Borel) version of QCD sum rules for fp [10]: 



fpm 4 



m 



Q 



where 



so 

exp (-a(mp - ttiq)) = J ds exp (-a(s - m 2 Q )) ^-Imll 5 \ Ap + n 5 (cr) | 



NP 



(10) 



U 5(cr)\ NP = -m Q (qq) I - 

| 8<J 7ra s p(qq) 
27 mQ 



+ 



l2imi Q (qq) ' 2 



1 - 



2 - 



2 6 



(11) 



Sum rules for the meson mass mp can be obtained from ratios of FESR with 
different n, or from the logarithmic derivative of (10) with respect to a. Appli- 
cations of other versions of QCD sum rules to the estimate of fp are presented 
e.g. in [11] and [2]. 

In (9) and (11), Mq parametrizes the dimension 5 quark-gluon conden- 
sate (g s qa ^G^ 'q) = M^qq), and p ^ 1 represents the deviation of the dimen- 
sion 6 four-quark condensate from pure factorization. Typical values, obtained 
from PCAC and from applications of the method to other channels, are e.g. 
{qq}(l GeV) ~ -0.016 GeV 3 , (a s G 2 ) ~ 0.02^0.06 GeV 4 , M 2 ~ 0.5^0.8 GeV 2 
and p ~ 1 -T- 3. In practice, results for fp are not sensitive to the value of the 
gluon condensate and p. Conversely, more important is their dependence on the 
input values of quark masses, typical values being mi,(Q 2 = m 2 ) ~ 4.6^-4.8 GeV 
and m c (Q 2 = m 2 ) ~ 1.3 -=- 1.4 GeV. 

Regarding the practical use of QCD sum rules, we remark that the non- 
perturbative expansions in (9) and (11) have been truncated to the first few 
significant terms (of lowest dimension) , neglecting higher condensates. The va- 
lidity of such an approximation clearly depends on the values of n or a. Indeed, 
larger n and a emphasize the ground state in (8) and (10) and minimize the 
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role of the continuum, but correspondingly also increase the size of higher power 
corrections. Thus, in practical applications one considers as a compromise the 
first few values of n, or a suitable range of a, such that the truncation of the 
OPE is justified. Another point concerns the role of the continuum threshold 
so, which is not determined a priori. Clearly, the predictions for the hadronic 
observables must not be sensitive to the values of this parameter. Thus, the 
procedure to exploit the FESR (8) is to search for a "duality window", i.e. a 
range in which the predicted tub (and mo) agrees with the measured one and 
is stable against s (and n). In this window one then solves for f B , and this 
estimate should be reasonably reliable. Similarly, in the case of the exponential 
sum rules (10) one imposes stability of the predicted mp and fp in the range of 
o and so mentioned above. 

In Table 1 we list the results for fg (and for completeness also for fp>) 
previously obtained from FESR and from exponential sum rules [12,13,14], using 
a s and (qq) evaluated at a scale uiq. The authors of [10] use nib = 4.8 GeV, 
which is somewhat larger than the values of m\, used by the other groups. This 
shows that the method is sensitive to the input value of the heavy quark mass. 
The "±" in Table 1 accounts for reasonable variations of the input parameters. 
Within SU(3) breaking effects, the results for fp> ~ fp> s are compatible with 
the experimental data [7]. 



Method 


/b//tt 


fo/U 


Ref. 


Laplace 


1.02 ±0.11 


1.33 ±0.19 


[10] 


Hilbert moments 


1.35 ±0.25 


1.70 ±0.20 


[12] 


Laplace & Hilbert 


1.38 ±0.14 


1.31 ±0.12 


[13] 


Laplace 


1.29 ±0.15 




[14] 


Laplace (HQET) 


1.48 ±0.39 


1.29 ±0.23 


[17] 


Laplace (HQET) 


1.7 ±0.2 




[19] 



Table 1 QCD sum rules predictions for fp 



Another important point is that the size of the 0(a s ) corrections in the 
sum rule for fg turn out to be rather large, of the order of 40%. Moreover, at 
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this two loop level there is an ambiguity in the definition of the argument of oi s , 
i.e. whether it should be or the much smaller confinement scale of the order 
of 1 GeV, which can make a substantial numerical difference. 

Indeed, in this regard it is instructive to consider the "static approximation" 
of (10), where tuq becomes infinitely heavy. This can formally be obtained by 
defining the "non-relativistic" variables: 

y/s = m Q +uj] m P = rriQ + uj r ; t = 2m Q a, (12) 

where r is a euclidean time and the "binding energy" ujr should be flavor- 
independent in the heavy quark limit. Expanding up to the first order 1/toq, 
and retaining the condensates up to dimension 5, Eq.(10) becomes [15,16,17]: 

(fpVmp~f ( — J exp {-tu r ) = 
\ m Q ) 







, . mn\ 13 47r 2 Auj 
2 In ( — — ) H H 

2u> / 3 9 tuq 



(qq) 



MqT / At \ : {a s G 2 ) 



1-— + 



16 \ tuq) l2imiQ(qq) 



(13) 



Analogously to the relativistic version (10), a sum rule for ujr can be obtained 
from the logarithmic derivative of Eq.(13) with respect to r. 

By construction, this non-relativistic form should reproduce the large jtiq 
dimensional behaviour implied by (2), plus 1/mg corrections. In addition, it 
should be noticed that, for tuq — > oo, the right hand side of (13) has no explicit 
reference to ttiq, and actually now involves mass scales of the order of the con- 
finement scale ~ 1 GeV (the same is expectedly true for ujr). This fact naively 
suggests to take a low energy, confinement scale as the argument of the running 
a s and (qq).* However, Eq.(13) clearly shows that there is a difficulty with the 
0(a s ) correction, associated with the appearance of a large (diverging) InmQ, 
so that the limit uiq — > oo is only defined after a resummation of these terms to 
all orders. Also, one can see that the finite coefficient in the a s correction is of 
the order of 10, giving rise to a large contribution, regardless of the choice of the 
argument of the strong coupling constant. The separation of the lnmg depen- 
dence from the sum rule and the definition of the appropriate scale require the 



* This suggestion was also made in [18], discussing the application of QCD sum 
rules to quarkonium states. 
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renormalization group improved summations of the leading (a s lnmg)" and the 
next-to-leading a s (a s lnmg)" terms to all orders in perturbation theory. This 
has been recently done by formulating the QCD sum rules in the framework of 
the HQET [17,19]. 

Basically, in this theory one has to investigate the two-point correlator 

U 5 (w) = i J d 4 xexp(ikx) < 0|T (j 5 (x) J 5 (0)^ |0 >, (14) 

where J5 = 975/1(3 (u) is the effective pseudoscalar current, with Hq(v) the 
velocity-dependent heavy quark field. The ground state contribution is essen- 
tially determined by (0| J§\P{v)) = v /mQi ? ( / u), where F(fj,) is an effective theory 
low energy constant, independent of uiq but depending on a renormalization 
scale fx. The relation between the physical and the HEQT fp is: 



f Pv ^ = C F (^)F(iu) + 0( A 



QCD 



(15) 



fi m Q 

where Cf is a QCD short-distance coefficient, which is computed up to the 
next-to-leading order [20,21,17]. The left side of (15) must be /x-independent; to 
the considered order the rriQ- and the /(/-dependences in Cf( — — ) factorize, so 

that (15) can be suitably rewritten as : 



f P y/mp~ = C F (m Q )F, 



(16) 



where 



C F (m Q ) = (a s (m Q )) 



d/2 



1 a s (m Q ) z 



n 



(<*s(fJ>))' 



-d/2 



TV 



F{H) 



(17a) 



(17.) 



In (17) d = 70/ A) with 70 and flo the leading order anomalous dimension and 
QCD /3-function, Z is a scheme-independent constant and 8 is scheme-dependent 
(for the explicit values we refer e.g. to [17]). According to (17&), F is a scale- 
independent, universal constant of QCD. 

By the same procedure leading to (13) one obtains a QCD sum rule of the 
form {rriQ = 00): 



F 2 (fx)exp(-ruj R ) 



LO 



du u) 2 exp (—tuj) < 1 + 



a 



TT 







4tt 2 

H — 2 In (tuj) + 5 
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M 2 r 2 
16 



21n( 



+ 



fJ,T , 



13 

y 



(18) 
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where the ellipses denote contributions of order a s and of higher dimensional 
condensates. Eq.(18) corresponds to (13) at tuq = oo, except from the replace- 
ment of vriQ by \i in the log and the appearance of 5. Both the /(/-dependence 
and 8 must cancel when considering the sum rule for the renormalization-group 
invariant F. Indeed, the summation of the ln/i to all orders produces the factor 
(a s (2/r)/a s ( / u)) d and, taking into account (17&), the running a s must be eval- 
uated at a low-energy scale of the effective theory. The structure of the QCD 
sum rule for F is finally [17]: 

(a s (2/r)f (l + 2^Z) F 2 exp (-ru R ) 

= ^ du u) 2 exp -tu H 1 2 In (too) 

TV 2 J 7T \ 3 9 / 



- (qq)U/t 



16 



(19) 



where a s is evaluated at a confinement scale. 

The analysis then proceeds by looking for ranges in r and in uq where the 
prediction for F is stable. The results exhibit a dependence on the value of lor, 
which is typically of the order of 0.5 0.8 GeV ', depending on the value of the 
heavy quark mass m& (or m c ). Such a dependence somehow limits the accuracy 
of the predictions. Also, there is an uncertainty associated to the actual choice 
for the low energy scale in the running strong coupling constant. Although the 
different choices formally have an impact at the a 2 level, in practice they make 
a substantial difference because the a s corrections are large by themselves. 

Combining the various uncertainties, the final result for the "static limit" 
of fs, defined as 

f stat _ C{m b )F 

is fg at = 200 -J- 300 MeV , in good agreement with the value recently found from 
lattice calculations [22]. However, the price is represented by a 100% order a s 
correction, which questions the validity of the perturbative expansion in this 
case. 

Concerning the 0(1 /toq) corrections due to the breaking of the heavy quark 
symmetry, a systematic analysis is possible by combining (19) with (13), and 
reintroducing the terms which vanish in the infinite mass limit [17,19]. One 
source of corrections, which can be immediately assessed, is related to the factor 
(mp/mo) on the left hand side of (13). This shows that, although naively of 
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order (Aqcd/itiq), the expansion parameter is actually a factor times (uJu/niQ). 
For current values of quark masses, 1/toq corrections are then expected to be 
reasonably small for fs, but very large for fjy. 

The results of the complete analysis of the sum rules including I/tuq cor- 
rections are reported in Table 1, and one can compare them with the values 
previously found from their relativistic counterparts. Within the uncertainties, 
the overall picture indicates a leptonic constant of the B comparable (or even 
larger) to that of the D, and in the range of possible measurements at a beauty 
factory. 

3. A determination of the beauty quark mass 

In the previous section we mentioned the dependence of the QCD sum 
rule estimates of fs on the value of m&. This is true also for other transition 
matrix elements. Thus, a precise determination of m& is highly desirable in 
order to reduce the uncertainty in the theoretical predictions. The ideal source 
of information on mj should be represented by the application of QCD sum 
rules to the T system, where the data are extremely precise. A typical result 
of such analysis, using exponential moments in the static approximation, is 
nib = 4.65 ± 0.05 GeV [23]. It turns out that estimates of nib are affected by 
the approximations used to confront theory and experiment, and in particular 
by the role of the a s corrections [2,23]. 

In [18] the convenience of studying ratios of the exponential moments (5) 
was pointed out. In fact, in ratios the importance of higher order radiative 
corrections could be reduced considerably, in particular in the non-relativistic 
limit. The latter corresponds to infinite quark mass, and only the leading mass 
term was considered in [18]. Moreover, since the ratios are found to depend 
sensitively on the quark mass, they should be suitable for the extraction of 
nib- As an attempt to pursue improvements over the past analyses, it should 
be interesting to reconsider the QCD sum rules proposed in [18], including the 
next-to-leading mass corrections to the non-relativistic limit [24]. 

To this aim, following the familiar procedure, we start from the correlator: 




(21) 



<"> = ;/ 

n Jo 



oc 



ds exp (—as)Imll(s), 



(22) 
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we define the ratio, sensitive to the mass of the ground state: 



BW = _>n W1 = 5V), 



(23) 



which tends to m\ for a — > oo. 



In the non-relativistic limit, ^fs = 2mb + u>, Eqs.(22) and (23) become 
respectively: 

1 f°° 

II(r) = — / du exp (— ru)lmll(uj), 
n Jo 

where a = r/4m&, and 



(24) 



R(a = r/4m b ) = 2m h - -^[lnU(r)]. 
At the two-loop level in perturbative QCD [25]: 



(25) 



-ImU(q 2 ) AF = -L v (3-v 2 ) 

TV (ill 



x 1 + 



4a, 



TV (v + 3) / 7T 3 

2v ~ 4 I 2 ~ 47r 



x 6(s - 4ra&), 



(26) 



where i> = a/1 — 4m^/s (m/, is here the 6-quark pole mass). The leading non- 
perturbative contribution is given by the gluon condensate [1]: 



U(q 2 ) NP 



48q 4 



3(v 2 + - v 2 ) 2 ^1 + v 3v 4 -2v 2 + 3 



2v 5 



1-v 



(^±G 2 ). (27) 

TV 



Expanding (26) to leading plus next-to leading orders in mj — > oo, we finally 
obtain the following QCD expression for R(r) defined in Eq.(25): 



R{t) =2mJl + 



1 / TV 



QmbT I 3 \rribT 



1/2 



a«(r) 



x 



2/19 3 5 

1 —TV H 

TvmbT \ 2 167T 8 



(28) 



where 



a s {r) = 



12tv 



(29) 



231n(4m b /rA 2 )' 

The ratio (28) must be compared with a corresponding ratio involving the ex- 
perimental data. We parametrize the latter by the sum of narrow T resonances 
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followed by a hadronic continuum starting at a threshold sq, modelled by per- 
turbative QCD as given in (26): 

U(a) E xp = j 2 — ^2 r v m v exp (-am v ) 

47r a EM v 

1 [°° 

H — / (is exp {—<js)ImTi(s)QCD- (30) 

Analogously to (13), one can notice from (28) that in practice the sum rule 
determines the difference between the hadronic ground state mass my and 2m;,, 
which is related to scales relevant to the confinement dynamics. 

The values of m& are then determined by matching the theoretical expres- 
sion (28) to the square root of the experimental ratio 

R{a) EXP = -^\nn{a) EXP . (31) 

Stability requires the matching to occur in a range of a and sq (actually the 
sum rule is extremely well saturated by the T family, so that the continuum 
has very little role), with corrections terms in (28) at a safe, few percent level. 
Varying Aqcd and of (a s G 2 ) in the currently allowed range one finds: nib = 
4.71 ± 0.05 GeV. 

As a comment to this determination, we notice that the subleading quark 
mass corrections can be important. For instance, the 0(a s /mb^Jmb) term in 
(28) is of the same size and sign as the non-perturbative term, for the relevant 
values of r. Hence, it would not be fully justified to keep the latter and ignore the 
former. Thus, at this level, with all correction terms safely small, the accuracy 
on m6 is essentially limited by the uncertainties in Aqcd, known to a factor of 
2, and in the gluon condensate, known to a factor of 3. 

4. Semileptonic B decays into charmed final states in the limit of 
infinitely heavy quarks 

The hadronic matrix elements that describe the decays into negative-parity 
states: 

B -> Dlv (32) 

and 

B -> D*tv (33) 
can be written in terms of form factors as follows: 

< D(p')\ c^b\ B(p) >= G+ ( P + P ')^ + G- </ M (34) 
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- < D*(p',e)\ B(p) >= i F e^ aP e*» (p + p')^ (35) 

- < D*(p', e)| c W b\ B(p) >= F e ; + F+ (e* .p)(p+p')^ + F. (e*-p)q^ (36) 

where G±, F, F ,± are functions of q 2 = (p — p') 2 . 

These form factors will be discussed for finite values of the heavy quark 
masses in Section 6. Here we consider the results obtained in the limit tuq — > oo. 
In the infinite heavy quark mass limit one can build up an effective field theory 
(HQET) [5] and obtain useful predictions for the above form factors. More 
precisely, one assumes that tuq is much larger than the average momentum of the 
light degrees of freedom q^: tuq » \q^\. If we write Pq = rriQV^ + k^ = rriQVQ, 
where is the hadron velocity, Vq is the quark velocity and /c M is a residual 
small momentum, in the limit 

niQ — > oo , Pq/tuq fixed (37) 

we get 

= (38) 

i.e. the heavy quark velocity is equal to the hadron velocity: QCD interac- 
tions do not change Vq that is therefore conserved. Eq.(38) is called velocity 
supers election rule [5]. 

On can show that in the limit (37) the heavy quark propagator becomes: 

t±li- (39) 
and the heavy quark gluon vertex assumes the form: 

-igT a v^ , (40) 
because, as a consequence of the equations of the motion, the effective heavy 

i> + 1 • 

quark field operator h v (x) = — - — e imQVX Q(x) satisfies the equation 

i)h v = h v . (41) 

These results follow from a strong interaction effective Lagrangian which is 
the starting point of HQET: 

£ = £^=]D ^D^K . (42) 
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The effective Lagrangian (42) has two symmetries that are not present in 
the original QCD lagrangian and are a consequence of the heavy quark mass 
limit: 

1) SU(2) spin symmetry, which follows from the absence of 7 matrices in (42). 

2) SU(N) (N = 2 if only charm and beauty are considered) flavour symmetry, 
which follows from the independence of the lagrangian from the heavy quark 
flavour. 

Similar to Eq.(2) relating fs to /d, because of these symmetries one can 
relate several physical quantities that are in principle independent. In particular, 
all the form factors in Eqs. (34-36) can be expressed in terms of one universal 
form factor £ (Isgur-Wise function). £ is a function of w = v ■ v', where v and v' 
are the B and D four velocities: 

£ = £H (43) 

with 

, m 2 B + m 2 D - q 2 . . 

w = v . v ' = — ° ±> — z_ . (44) 

2mBTnD 

Notice that £(1) = 1 since the vector charge is the generator of the flavour 
symmetry. 

The matrix elements in (34-36) are expressed in terms of £ (neglecting 1/toq 
corrections) as follows: 

< D{v')\V»\B{v) >= ^m B m D £(w) (v + v'Y (45) 



< D*(v',e)\V - A^\B(v) > = ^m B m D £(w 

- (1 + w)e^ + (e* • v)v'" 



■ fj.va/3 * 1 
It t v U a Up 



(46) 



The relations among the form factors in Eqs. (34-36) and the universal func- 
tion £ can be obtained by a straightforward calculation. 
We shall also consider the semileptonic transitions 

B -> D**tv , (47) 

where D** are positive parity p-wave charmed mesons. According to the quark 
model one expects four states: 2S+1 Lj = 3 P 2 , 3 -Pi, 3 -Po (corresponding to total 
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spin S = 1) and 1 P± (corresponding to total spin S = 0). In the HQET the 
heavy quark spin sq is decoupled from the light quark angular momentum sg = 
J — sq, so that it is convenient to work with the two degenerate multiplets 
J p = (0^y 2 , l^/ 2 ) an d J P = (I3/2' ^3/2) wn i cn differ by the value sg = 1/2 
and 3/2 respectively (the relative orbital momentum is L = 1). The 1 + states 
defined above are given by the following linear combinations [26,28]: 

\^ / 2>=^l\ 1 Pl>+Jl\ 3 Pl> 



(48) 

while \0y 2 >= l 3 -^" > an< ^ 1^3/2 >= l 3 -^" 2 >- ^he 2y 2 state has been observed 
experimentally with a mass of 2460 MeV: it is denoted by £2 (2460). The l+ /2 
meson can be identified with D!(2420) even though a ly 2 component is likely 
to be contained in this state [26,29]. They are both narrow (r < 20 MeV) since 
their strong decays occur in d-wave, in contrast with the states D$ (the 0^ 2 
state) and D[ (the 1^ 2 state) that can also decay in s-wave. 

In the HQET the matrix elements for the decays (47) can be written in 
terms of two universal functions Ti/ 2 (w) and t 3 / 2 (w) [27]. The form factor 
T3/ 2 (w) describes the transitions into the 2^ 2 and 1^ 2 states, and Ti/ 2 (w) those 
to the 1^2 an d 0]y 2 states. The hadronic matrix elements for the decays (47) 
in the tuq — > 00 limit can be written as [27]: 

< Doiy'^A^Biy) >= -2^m B m D ** r 1/2 (w) (v - v% (49) 

< D[(v', e)|V^ - A^\B(v) >= y/m B m D .* t 1/2 (w)x 

x [2 (w - 1) e* + 2 (e* • 1;)^ - ze^ a/3 e*> + </)> - u')"] 

(50) 

< D x {v\ e)\V^ - A^\B(v) >= y/m B m D ** r 3/2 (w)x 

< D 2 (v', e)|V M - A^\B(v) >= y/m B m D ** t 3/2 (w) 

[ i ^e^e^v^v + v'f(v - v'f - V3(w + l)e> Q + (52) 
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where now = c v i^^b v , = <V7 M 756«, c v > and b v are charm and beauty 
quark operators in HQET, and m^** is the mass of the positive parity charmed 
resonances. In the following Section we shall review the QCD sum rules calcu- 
lation of £(w), Ti/ 2 (w) and t 3 / 2 (w). 

5. Three point functions 

To obtain QCD sum rules for the form factors one starts from the three- 
point correlators [30]: 

H r ^(p,p', q ) = i 2 J dxdye'V-*-™) < 0\T {j r (x)J^ A (0)J 5 (y)) \0 > . (53) 

In (53) Jr denotes the currents interpolating the final charmed meson with T 
the appropriate set of Lorentz indices (e.g. Jr = qc, qi^c and q(^ u d ^ +7 M d „ — 

\9nvl p d p )c for ± and 2+ respectively), and J5 is the analogous operator for 
the initial B- meson, see Eq.(3). The form factors are extracted by expanding, 
in each case, Eq.(53) in independent Lorentz structures. 

Similar to the analysis of two-point functions, we write a sum rule by com- 
puting (53) in two different ways: by including hadronic states plus a continuum 
of states modelled by perturbative QCD on one side, and by performing an OPE 
in the framework of QCD on the other side, accounting for the asymptotic free- 
dom contribution plus higher dimensional non-perturbative condensates. Fi- 
nally, these alternative representations are matched to one another, and the 
convergence is improved by one of the methods discussed in Sect. 2, e.g. by 
taking the Borel transform of both sides of the sum rule. The double Borel 
transform, as applied to any of the Lorentz invariant functions Ui(p 2 , p' 2 , q 2 ) 
appearing on the left hand side of (53), looks as follows: 



Ui(a, a', q 2 ) = J Pi{ s -> s ' \q 2 ) exp(-sa) exp(-sV) ds ds' , (54) 

where IL(cr, a', q 2 ) is obtained applying the Borel transform to the OPE of 
Ui(p 2 ,p' 2 , q 2 ), and the double spectral function pi(s,s',q 2 ) is saturated by the 
lowest hadronic states, involving the form factors of interest, plus the QCD 
continuum starting at some thresholds so, s' . 

In order to perform the heavy-quark mass limit and then estimate the de- 
pendence of the the form factors on mi, we estimate the sum rules at q 2 = 0, 
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and take the limit tuq = mi, — > oo with m c = rrib/vZ, Z fixed [31]. The form 
factors then become functions of w via Eq.(44) at q 2 = 0, which reads: 



\[Z = w + \J w 2 — 1 



(55) 



In the case of 0~ — > (0~, 1~) semileptonic transitions, this procedure reproduces 
the main features of the results obtained by directly working in the infinite 
heavy-quark limit [17,32]. As for the Borel parameters a and a', in the heavy- 
quark limit, in analogy with (12), we write T\ = 2 a and = 2 m c a'. For 
the continuum thresholds we put s' = m 2 /Z + 2 [i+/y/~Z, so = ml + 2 m& \i- 
and for the meson masses we introduce the "binding energies" as follows: m 2 B — 
m 2 b + 2m b ur, m 2 P = m 2 h /Z + 2 m b u R ' '/y/Z where P = D,D*, D**. 
The resulting sum rules for the universal form factors are [33] : 



-<qq> 



r 3 Z 

^14^2 (v^+l)3 

A 2 r 2 («;-1) 



H(w) 



1 - 



48 



+ 1)/16)| 



(56) 



Tl/2(w) = 



f(+)F I- 8^2^+1 



J( W ) 



< qq > 



A 2 TiT 2 

12 



(w + l) 



e -(A 2 (r 1 2 /8+r 2 2 /8+2»r 1 r 2 /8))| ^ 



e (^HT 1+ a;^T 2 ) ^ Z 3/2 
T3/2{W) = — — ; j= ; — J(w) , 



(58) 



/ T F 7T 2 (v^+l)5" 

where A 2 = Mq/2 arises from the contribution of the dimension 5 operator 
introduced with regard to Eqs.(9) and (11). In Eq.(56) we have set T\ = T2 = r/2 
and ur = u>' R , which ensures £(1) = 1 for any value of r. Actually, once 
F is replaced by the sum rule (19), the dependence on ujr drops from the 
ratio, and this minimizes the sensitivity of £(w) on the values of this quantity. 
Explicit expressions for 0(a s ) corrections to (56) have been presented in [34]. 
The complete a s corrections for the case of (57) and (58) are not available yet. 
In (57), (58) we use the leptonic constants and fr defined as follows: 



<0|^|M(0+ ;p) >=i- 



+ 



L/2' — > 

< 0|J^|M(2+ /2 ;p,e) >= e^fry/m^ 



(59) 
(60) 
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The numerical values for fr and are obtained, in analogy with F, by 
two-point QCD sum rules [33] and the results are as follows: fa = (0.18 ± 
0.02) GeVi and /<+) = (0.46 ± 0.06) GeVi . 

The functions H(w), J(w) and I(w) in (56)-(58) are given by: 



#(«,) = — L-/' " dy [ ~ dy' (y + y>) e-^y+y'^ 4 0(2 wyy'-y 2 -y 
V Z — 1 Jo Jo 



(61) 



J(u,) = (v ^_ 1)3 J o dy jf dy' (y - y^^y x 

O^yy'-y 2 -?/ 2 ) (62) 



7H = — J- ^ ^'[(y 2 + 2yy0 (Z - Vz + 1) - Sv^y' 2 ] x 



e -( y r 2 /2+ y 'r 1 /2) e(2wyy /_ 2/ 2_ 2/ /2 ) ^ (63) 



where are the effective thresholds in the heavy quark limit. 

The stability analysis of the results for the form factors is an obvious ex- 
tension of the method described in Sect. 2. One observes explicitly that the non- 
perturbative terms turn out to be numerically negligible for r 3 / 2 , and therefore 
have not been included in (58). Moreover, it should be noted that the optimal 
values of the Borel parameters T\ and r 2 as well as the continuum thresholds fx + 
and fj,- are taken at the same values obtained in the analysis of two-point func- 
tions [33], in particular \i_ = 0.7 ± 0.1 GeV and (1+ = 1.3 ± 0.1 GeV. Finally, 
the form factors T\/2 and T3/2 should include a multiplicative 0(a s ) correction 
which depends on the anomalous dimension of the current operators and on w. 
Its numerical effect is small and can be included in the theoretical uncertainties 
of the results. Furthermore, as far as £ is concerned, we notice that the large 
0(a s ) corrections affecting F (see discussion in Section 2) should drop from the 
ratio in (56), which also preserves £(1) = 1. 

An important role in the numerical analysis of £ is played by the continuum. 
In spite of the denominators (\/~Z — l) n = (w — 1 + Vw 2 — l) n , with n = 1, 3, the 
integrals H(w), J(w) and I(w) are finite for w — > 1. On the other hand the slope 
of the universal function £ turns out to be divergent, a result which is unnatural. 
As shown in [17] this is an artifact of the approximations used in modelling the 
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higher mass resonances as a continuum of states. To cure this problem for 
the "Isgur-Wise" function £(w) one can modify the integration region in (61) 
without affecting the normalization constraint at zero recoil £(1) = 1 (which is 
easily seen to be verified by the sum rule). One can replace, for example, the 
6 function by 6(2 wyy' — y 2 — y' 2 ) x 6[2/i(l + w — \Jw 2 — 1 — y — y'\ (case 
a) or, as another possibility, one can just change the upper limit of integration 
2/i — > 2/x(l + w + y/w 2 — + w) (case b). Both procedures lead to a finite 
slope. However, they add some significant uncertainty to the determination of 
£ (w) for values of w relevant to nonzero recoil, and in particular by varying the 
continuum integration domain in Eq.(56) one obtains different values for the 
slope £'(1). The value that better fits the experimental data corresponds to the 
case a and is [17]: 

£'(1) = -1.28 ±0.25 , (64) 
while in [35] it is suggested to adopt the case b, which leads to: 

= -0.70 ±0.10 . (65) 

An analogous result has been found in [34] including 0(a s ) corrections. We 
observe that all these values satisfy the lower bound: 

|£'(1)|> 0.25 (66) 

obtained in [36]. 

Regarding the numerical evaluation of Ti/ 2 (w) and r 3 / 2 (w) according to 
Eqs.(57)-(58), we notice that, in this case, there is no general normalization 
condition on the form factors at w = 1. In Fig.l we depict the functions T\/2{w) 
and t 3 / 2 (w) for the central values of // + and \i- and with 1/ri ~ I/V2 — 1.5 GeV. 
We observe that the numerical values at w = 1: Tl/2(1) — 73/2(1) — 0.22 are in 
reasonable agreement with the outcome of the constituent quark model in [27]. 
Moreover, they also agree with the estimate 7~i/ 2 (l) — 0.35 ± 0.2 given in [35]. 
Finally, we observe that the values in Fig.l are inside the upper bound put by 
the optical sum rule [37]: 

|r 1/2 (l)| 2 + 2|t 3 / 2 (1)| 2 < \ m °~ mc c 0.5 . (67) 

In Table 2 we list the values of the branching ratios resulting from the 
analysis presented in this Section. 
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Channel Branching Ratio 



B^D(l-)£v 4.3 xlO" 2 (^) 

B^D £u 5xl0" 4 (^) 2 

B^Dltv 7xl0" 4 (^) 2 

B^Dytv 1 x 10" 3 (^) 

B ->D**£v 2 x 10" 3 (^) 2 



Table 2 . Branching ratios obtained for the various B-meson semileptonic transitions 
in the heavy-quark limit [33] (tb = 1.21 ps). 

6. Semileptonic form factors for finite heavy quark masses 

The method of QCD Sum Rules can be applied to the evaluation of the form 
factors in Eqs. (34-36) without performing the limit uiq — > oo. This offers the 
possibility to evaluate systematically 0(l/mg) corrections for the transitions 
B — > D, D*. We shall consider also the form factors for decays into the positive 
parity 0^ 2 and 1^ 2 states (the other form factors have not been computed yet). 
The form factors are given by: 

- < Dq*(p') \ c 7m75 6| B(p) >= g+ (p' +p )^ + g_ % (68) 
- < D[(p',e)\ c 7m 75^| B(p) >= i f e^e'V +P)V (69) 

- <D' 1 (p',e)\c^b\B(p)>= f oe ; + f + (e*-p)(p' + p)^ + f-(e*-p)q^ (70) 

where e M is the D[ polarization vector and the form factors g±, f, /o, f± are 
functions of q 2 . For this application it is also necessary to evaluate via QCD 
sum rules the vacuum-to-meson constants: 

< 0\qc\D o >= f .2^<L (71) 
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and 

< 0\q ltxl5 c\D[(p, e) >= . (72) 

9d[ 

The result for the mass is m£> = rn^ = 2.5 ± 0.1 GeV , and the leptonic 
constants are given in Table 3 [38]. 



Leptonic constant 


Theoretical result 


fD 


170 ± 20 MeV 


9d[ 


9.8 ±1.5 



Table 3 Leptonic constants, see Eqs. (71), (72). 



Our results for the values of the various form factors at Q 2 = are reported 
in Table 4. As for the <5 2 -dependence of the form factors, it could be derived 
by QCD sum rules; however, for simplicity, we assume a simple pole behaviour 

Q 2 

Fi(Q 2 ) = Fj(0)/(1 H j), which is adequate in the present case. This is also 

suggested by dispersion relations, where the pole masses are those of be mesons 
with appropriate quantum numbers. 







9+(0) 


/(0) 


/+(0) 


/o(0) 






G+(0) 


F(0) 


F+(0) 


^o(O) 








(GeV -1 ) 


(GeV -1 ) 


(GeV) 


B - 


-D(O-) 


0.67 ±0.19 








B - 


■+ D*(l~) 




0.12 ±0.02 


-0.05 ±0.02 


4.35 ±0.15 


B - 


-> D^(0+) 


0.30 ±0.03 








B - 


■+ D**(l+) 




0.10 ±0.03 


-0.02 ±0.02 


0.9 ±0.3 



Table 4 Form factors for B decays for finite tuq 
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The results for the widths are rather insensitive to the precise value of 
these masses; we take m = 6.34 GeV for the 1~ pole (for G±, F, fo and / + ) and 
m = 6.73 GeV for the 1 + pole (for g + , /, Fo and F + ). 

From the values in Table 4 we can compute the branching ratios for the 
decays B — > D,D* , B — > Dq, D[. The results are collected in Table 5. 



Decay mode BR (Theory ) BR (Exp. [39]) 



B - 


-> Div 


1.5(fe) 2 xl0- 2 


1.75 ±0.42 ±0.35 x 10" 2 


B - 


-> DHv 


4.6(^) 2 xl0" 2 


4.8 ±0.4 ±0.7 x 10" 2 


B - 


■* D Q £u 


0.15(^) 2 xl0" 2 




B - 


-> D[£u 


0.15(^) 2 xl0" 2 





Table 5 Branching ratios for finite itiq 



Comparing these results with those in Table 2, we observe that the 1/mg 
corrections affect much more strongly the decays into positive parity charmed 
states. This can be qualitatively understood by recalling that in the case of 
higher charmed resonances the expansion parameter can be expected to be oc 
(mo** —rn c )/m c , which is a number of order one. For the transitions into D, D* 
the corrections linear in 1/mq are forbidden in the most relevant form factors 
at w = 1 [40] and this seems to protect the leading order result. Actually, the 
systematic analysis of B decays into negative parity states at the next-to-leading 
order in I/tuq implies the introduction of four new functions, in addition to £ 
[40] . Some of these functions are also subject to normalization conditions at the 
zero recoil point w = 1. They have been recently analyzed in the framework of 
QCD sum rules in [41]. 

7. Concluding remarks 

Summarizing the results presented in the previous sections, QCD sum rules 
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provide a conceptually consistent (and numerically simple) method to estimate 
from first principles semileptonic and purely leptonic transitions of the B meson. 

For the case of predicted values are in the potential of a B factory, and 
indeed such a measurement would be really welcome for theoretical reasons. 

For B — > D,D* semileptonic transitions, the method can provide a sys- 
tematic analysis of the relevant form factors, both for finite and for infinite m&. 
In the last regard, the results (and the simplicity) of the HQET can be easily 
incorporated in this framework, and the Isgur-Wise function £(w) can be pre- 
dicted, naturally accomodating the zero-recoil normalization £(1) = 1. Order a s 
effects as well as the form factors relevant to the general next-to-leading 1/toq 
structure can be also accounted for. Thus the situation with the transitions to 
negative-parity charmed states looks well settled from the theoretical point of 
view, waiting for more accurate experimental data, which would allow a precise 
and model independent determination of V c \,. 

Also the B — > D** semileptonic transitions can be evaluated in the frame- 
work of QCD sum rules, where D** represents the orbitally excited charmed 
resonances. Indeed, the analysis of these decays is theoretically related to the 
previous ones in this framework, which therefore provides a unified treatment 
of all form factors in terms of the same fundamental parameters. Moreover, the 
case of D** final states can be discussed in the view of the heavy quark mass 
limit and the two universal functions Ti/ 2 {w) and r 3 / 2 (w) can be predicted. The 
current analysis, more difficult as it involves operators with derivatives of the 
quark fields, is not complete as it is in the case of £(«;), in particular the order 
a s effects are not taken into account yet. In addition, the r's are not normalized 
at zero recoil so that, in general, uncertainties are larger and predictions can be 
considered as more qualitative. Nevertheless, it is encouraging that the values of 
Tj(l) thus found are in reasonable agreement with the independent quark model 
estimates, and comply with the theoretical upper bound from the optical sum 
rule. 

Concerning the branching ratio of B — > £)**, the results indicate values of 
the order of a few x 10 -3 for the sum of the various p-wave resonances, which 
is not enough to fill the observed gap between the inclusive semileptonic rate 
and the sum of the exclusive rates B — > D, D*£u. Clearly, this issue deserves 
further theoretical attention and represents an interesting point for the physics 
at B factories. 
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